Vector solitons in (2+1) dimensions 
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We address the problem of existence and stability of vector spatial solitons formed by two 
incoherently interacting optical beams in bulk Kerr and saturable media. We identify families 
of (2+l)-dimensional two-mode self-trapped beams, with and without a topological charge, and 
describe their properties analytically and numerically. 
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Recent experimental observations of multidimensional 
spatial optical solitons in different types of nonlinear 
materials [[[) call for a systematic analysis of the self- 
trapping of light in higher dimensions. When two (or 
more) fields interact nonlinearly, they can form multi- 
component trapped states, known as vector solitons. Vec- 
tor solitons, first theoretically studied in a (1+1)-D mod- 
els H|, were observed in birefringent fibers and planar 
waveguides fjj. Fabricated waveguiding structure local- 
izes such solitons in one of the two directions transverse 
to the direction of propagation, hence these solitons are 
effectively one-dimensional. It is only recently, that the 
theory and experiments on incoherent interaction and 
truly two-dimensional self-trapping of beams in a bulk 
(saturable) medium merged Q, indicating progress to- 
wards the observation and study of different types of 
(2+l)-D vector solitons and their interactions. 

The practical possibility of such an observation greatly 
depends upon the soliton stability in the media with re- 
alistic, Kerr or saturable nonlinearity. It is known that 
scalar (one-component), fundamental (2+l)-D solitons 
are stable in saturable media JlJ , but they exhibit critical 
collapse in Kerr-type media Jpf. However, as in the case 
of (1+1 )-D vector solitons |5[, both existence and sta- 
bility of multi-dimensional vector solitons are nontrivial 
issues, which have not been systematically adressed so 
far. 

In this Letter, we study (2+l)-D vector solitons in Kerr 
and saturable media. We analyze two classes of such soli- 
tons. First, we consider solitons formed by the coupling 
of two fundamental modes; such solitons are always bell- 
shaped. Secondly, we analyze the coupling between the 
fundamental mode of one field and the first-order mode 
(i.e. that carrying a topological charge) of the other field. 
In the latter case the vector solitons may possess a ring 
structure and are expected to be analogous to the two- 
hump (1+1)-D vector solitons recently proved to be sta- 
ble in a saturable medium jf|. 

We consider two incoherently interacting beams prop- 
agating along the direction z in a bulk, weakly nonlin- 
ear optical medium. For a Kerr medium, the problem is 
described by the normalized, coupled equations for the 



slowly varying beam envelopes, E\ and Ei, 
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where Aj_ is the transverse Laplacian, and a measures 
the relative strength of cross- and self-phase modulation 
effects. Depending on the polarization of the beams, the 
nature of nonlinearity, and anisotropy of the material, a 
varies over a wide range. For a Kerr-type material with 
nonresonant electronic nonlinearity a > 2/3, whereas for 
a nonlinearity due to molecular orientation a < 7 Q . 
We look for solutions of Eqs. (fi\) in the form 
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where /3i and /3 2 are two independent propagation con- 
stants, and mi. 2 = 0, ±1 are topological charges. Mea- 
suring the radial coordinate in the units of and in- 
troducing the ratio of the propagation constants, A = 
P2/P1, from Eqs. (||) we derive a system of station- 
ary equations for the radially symmetric, normalized en- 
velopes u and v: 



A r u — u + (u 2 + av 2 )u = 0, 
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V — Xv + (v 2 + OV?)v = 0, 



(3) 



where A r = (l/r)(d/dr)(rd/dr), and we assume mi = 0. 
Following the notations introduced in |J , we describe all 
vector solitons (^|) by their "state vectors" |m 1 ,m 2 ). 

First, we consider solutions |0, 0). The families of these 
radially symmetric, two-component vector solitons are 
characterized by a single parameter A, and at any fixed 
value of a, their existence domain is confined between 
two cut-off values, Ai and A 2 . When A < Ai or A > A 2 , 
self-trapping of coupled fields does not occur, and there 
exist only scalar solitons for either the u or v components. 
However, for Ai < A < A 2 , a two-mode self-trapped state 
emerges. Near the cutoff points (A w Ai or A w A 2 ) this 
state can be presented as a waveguide created by one 
field component with a small-amplitude guided mode of 
the other field component. Examples of |0, 0) solitons 
are presented in Figs. l(a, b) for a — 2. On the param- 
eter plane (cr, A) , the existence domain is skirted by two 
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curves, cti(A) and (72(A), which are defined by the corre- 
sponding cut-off values, Ai and A2, of the soliton-induced 
waveguides (see Fig. |l|). 
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FIG. 1. Existence region for the |0, 0} vector solitons. Solid 

- numerically obtained cutoff curves, o"i(A) and 02(A). Dashed 

- results of the variational analysis. (a,b) Amplitudes of the u 
(solid) and v (dashed) components of vector solitons at a — 2, 
corresponding to the points a and 6 on the plot above. 



For a = 1, the |0, 0) vector solitons exist only at 
A = 1, and their properties resemble those of the (1+1)- 
D Manakov vector solitons. They can be constructed 
by the transformation u = U cos 6 and v — U sin 9, 
where 9 is arbitrary and U satisfies the scalar equation 
d 2 U/dr 2 + (l/r)(dU/dr) - U + U 3 = 0. 

To describe the existence domain of the multi- 
dimensional vector solitons analytically, we employ the 
variational technique ||. We look for stationary two- 
component solutions of Eqs. (m in the form u{r) — 
Aexp{— r 2 /a 2 ), v(r) = Bexp(— r /b 2 ), where the param- 
eters A, B, a, and b are defined by variation of the ef- 
fective Lagrangian of the model (|l|) . Details of this anal- 
ysis will be published elsewhere. Here, we mention that 
the coupled algebraic equations of the variational analy- 
sis allow us to find the borders of the existence domain 
for the 1 0, 0) vector solitons: cti(A) = (1 + ^fX) 2 /A and 
0-2 (A) = (1 + \/A) 2 /(4A), shown in Fig. | by dashed 
curves. One can see that the variational approach pro- 
vides an excellent alternative to the numerics in identi- 
fying the existence domains of the vector solitons. 

An important physical characteristic of vector solitons 
of this type is the total power defined as P — P u + P v = 
2ir L (u 2 + v 2 )rdr, where the partial powers P u and P v 
are the integrals of motion for the model ([!]) . Figures 2 
(a,b) show the total power of the (2+l)-D vector solitons 
vs A for a = 2/3 and a = 2, respectively. 
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FIG. 2. Total power of the vector soliton (solid), and the 
partial powers of its components (dashed and dot-dashed). 

As follows from these results, the cases a > 1 and 
<T < 1 are qualitatively different. In the former case, the 
total power of the vector soliton is lower than the power 
of a scalar soliton, P ~ll-7atcr = A = l [see Fig. 
2(b)]. This is an important and unexpected physical re- 
sult indicating, in contrast to the commonly held belief, 
that the excitation of vector solitary waves would require 
lower input power in comparison to scalar solitons. For 
cr < 1, the situation is opposite, and the vector solitons 
exist at higher power than scalar solitons [see Fig. 2(a)]. 
A lower total power of the vector soliton, as compared to 
the corresponding power of a scalar soliton, allows us to 
explain an effective suppression of the blow-up instabil- 
ity numerically observed at a = 2, for the special case of 
A = 1 , when an analytical form of the vector soliton can 
be found by a Hartree-type ansatz ^ . 

Similar to the well studied case of the (1+1 )-D vector 
solitons (see, e.g., Ref. 0), the soliton-induced waveg- 
uide can guide higher-order modes. In higher dimensions, 
such higher-order modes carry a topological charge, i.e. 
ni2 7^ in Eq. (^). In the simplest case, we consider 
rri2 = ±1 and analyze the system of stationary equations 
for the radially symmetric wave envelopes, Eqs. (|3|). Ex- 
amples of two-mode |0, 1) solitary waves are shown in 
Figs. 3 (a,b). 
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FIG. 3. Intensity profiles of vector solitons (solid), formed 
by scalar components with (v, dash-dotted) and without (u, 
dashed) a topological charge (cr = 3). 

Near the cut-off of the first-order mode, the second 
component (v) appears as a guided mode (dash-dotted) 
of the effective waveguide created by the scalar soliton in 
the w-component (dashed) . The total intensity, shown in 
Figs. 3 (a) by a solid curve, has a maximum at the beam 
center. Akin to the humps in the total intensity profile 
of (1+1)-D solitons ||, the ring shape of (2+l)-D vector 
solitons develops far from the cut-off for the first-order 
mode, when the guided mode deforms the soliton waveg- 
uide and creates a coupled state that has the maximum 
shifted from the beam center [see Fig. 3(b)]. 

Next, we conduct the numerical stability analysis, and 
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an accurate check of the Vakhitov-Kolokolov stability cri- 
terion for vector solitons ||. Our analysis reveals that 
neither the fundamental nor the first-order mode of the 
soliton-induced waveguide can arrest the collapse of the 
scalar bright (2+l)-D solitons in a Kerr medium, and the 
corresponding vector solitons |0, 0) and |0, 1) are linearly 
unstable. 
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FIG. 4. (a) Existence domain for |0, 1) solitons in a sat- 
urable medium, (b) Evolution of the total intensity at r — 
for different values of s (A = 0.6 and z = 50). 
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FIG. 5. Typical evolution of the |0, 1) solitons in a sat- 
urable medium (s = 0.65, A = 0.6). Shown are (a) the inten- 
sity distribution at z = 0, (b) evolution of the intensity profile 
at y — 10, and (c) the intensity distribution at z = 100. 

Given the established stability of (1+1)-D vector soli- 
tons in a saturable medium H , a tempting task is to look 
for the |0, 1) ring-like solitons in such a medium. For the 
(completely solvable) model of the so-called threshold- 
type nonlinearity, |0, 1) vector solitons have been recently 
analyzed in |l(| . Here we consider the model correspond- 
ing, in the isotropic approximation, to the physically 
realised solitons in photorefractive materials. The nor- 
malized dynamical equations for the envelopes of two in- 
coherently interacting beams can, in this case, be ap- 



proximately written in the form: 



' dz 
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Eia (1 + |^i,2 1 2 + |^2,i I 2 ) = 0- Seeking stationary so- 
lutions in the general form (Q) , and introducing the rela- 
tive propagation constant A = (1— 02)10-— Pi), we arrive, 
after corresponding renormalizations || , to the following 
system of equations [cf. Eqs. (g)]: 

A T u-u + uf(I)=0, (4) 



A r v ±v 



Xv + vf(I)=0, 



where /(/) = 1(1 + si)' 1 , I = u 2 + v 2 , and s = 1 - Pi 
plays the role of a saturation parameter. For s = 0, 
this system describes the Kerr nonlinearity (with a = 1). 
Since the contributions from the self- and cross-phase 
modulation are, in this case, equal, the lowest-order bell- 
shaped |0,0) solutions only exist at A = 1. In the re- 



maining region of the parameter plane (s,A), the solu- 
tions |0, 1), similar to those described above for the Kerr 
nonlinearity, are found [see Fig. 4(a)]. Again, the ring- 
shaped structure of these solutions develops only far from 
the cut-off for the vortex-type guided mode. Close to the 
cutoff, all vector solitons are bell-shaped. 

Although our numerical simulations confirmed that the 
saturation does have a strong stabilizing effect on the 
|0, 1) vector solitons [cf. cases s = and s — 0.6 in Fig. 
4(b)], vector solitons of this type appear to be linearly 
unstable. The instability, although largely suppressed by 
saturation, triggers the decay of the solitons into a dipolc 
structure (as shown in Fig. 5) for even small contribution 
of the charged mode. However, the quasi-stable dynamics 
exhibited by these vector solitons over long propagation 
distances may have serious implications on the attempts 
to observe them experimentally. Indeed, since, for the 
current experiments on spatial solitons in photorefractive 
media 11 1, the propagation length z = 100 corresponds 
to a crystal length of up to 40 mm, this slow-developing 
dynamical instability may be hard to definitively detect 
in experiments. 

In conclusion, we have analyzed multi-dimensional vec- 
tor solitons composed of two incoherently coupled beams 
in both Kerr and saturable media and defined, analyti- 
cally and numerically, existence domains for new classes 
of spatial solitary waves. A stabilizing effect of nonlinear- 
ity saturation on (2+l)-D vector solitons with topological 
charge has been demonstrated. 

Yu. S. K. thanks M. Segev for useful discussions and 
a copy of Ref. [p0| prior to its publication. 
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